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Abstract 

Nonstationary Gabor frames were recently introduced in adaptive signal analysis. They 
represent a natural generalization of classical Gabor frames by allowing for adaptivity of 
windows and lattice in either time or frequency. In this paper we show a general existence 
result for this family of frames. We then give a perturbation result for nonstationary Gabor 
frames and construct nonstationary Gabor frames with non-compactly supported windows 
from a related painless nonorthogonal expansion. Finally, the theoretical results are illus- 
trated by two examples of practical relevance. 

Keywords: adaptive representations, nonorthogonal expansions, irregular Gabor frames, 
existence 



1. Introduction 

The principal idea of Gabor frames was introduced in [H] with the aim to represent signals 
in a time-frequency localized manner. Since the work of Gabor himself, a lot of research has 
been done on the topic of atomic time-frequency representation. While it turned out that 
the original model proposed by Gabor does not yield stable representations in the sense of 
frames [HEHTUj, the existence of Gabor frames was first established in the so called painless 
case, [7j, which requires the use of compactly supported analysis windows. The existence of 
Gabor frames in more general situations was proved later [191 E3] and the proof often uses 
an argument invoking the invertibility of diagonally dominant matrices. 

Various irregular and adaptive versions of Gabor frames have been introduced over the 
years, cf. [TJ El [T^J |2T]. In these approaches, the irregularity usually concerns either the 
sampling set, which is allowed to deviate from a lattice, or the window, which is allowed 
to be modified. In [lj, varying windows as well as irregular sampling points are allowed, 
however, existence of a local frame is assumed, from which a global frame is constructed. 
Nonstationary Gabor frames give up the strict regularity of the classical Gabor setting, but, 
as opposed to irregular frames, maintain enough structure to guarantee efficient implementa- 
tion and, possibly approximate, efficient reconstruction. In analogy to the classical, regular 
case [7], painless nonstationary Gabor frames were introduced in [2], where the principal 
idea is described and illustrated in detail. The construction of painless nonstationary Gabor 
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frames is similar to, but more flexible than the construction of windowed modified cosine 
transforms and other lapped transforms [17] EI] that allow for adaptivity of the window 
length. An efficient and perfectly invertible constant-Q transform was recently introduced 
using nonstationary Gabor transforms [22] • In this and similar situations, redundancy of 
the transform is crucial, since non-redundant versions of the constant-Q transform lead to 
dyadic wavelet transforms, which are often inappropriate for audio signal processing. 
Redundancy allows for good localization of both the analysis and synthesis windows, and 
their respective Fourier transforms and often promote sparse representations in adaptive 
processing. 

Painless non-orthogonal expansions can only be devised if the involved analysis windows 
are either compactly supported or band limited. This requirement may sometimes be too re- 
strictive. For instance, one may be interested in designing frequency-adaptive nonstationary 
Gabor frames with windows that are compactly supported in time, i.e. can be implemented 
as FIR filters, cf. [11] and lend themselves to real-time implementation, cp. [9]. 

The present contribution addresses the case of nonstationary Gabor frames with more 



general windows than used in the painless case. In Theorem 3.4, we derive the existence 



of nonstationary Gabor frames directly from a generalized Walnut representation: under 
mild uniform decay conditions on all windows involved, we show an existence result of 
nonstationary Gabor frames in parallel to the result given in [23] for regular Gabor frames, 
also compare [T51 Theorem 6.5.1]. Note that the existence of a different class of nonstationary 
Gabor frames, the quilted Gabor frames, [8], was recently proved in the general context of 
spline type spaces in the remarkable paper [T8] . 

This paper is organized as follows. In the next section, we introduce notation and state 
some auxiliary results. In Section |3j we first define nonstationary Gabor frames and recall 
known results for the painless case. In Section 3^ a Walnut representation and a correspond- 
ing bound of the frame operator in the general setting is derived and Section |3~3] provides the 
existence of nonstationary Gabor frames. In Section 3^ we pursue two basic approaches for 
the construction of nonstationary Gabor frames . Using tools from the theory of perturbation 



of frames, we construct nonstationary frames from an existing frame in Proposition 3.7 In 



Corollary 3.8 we design nonstationary Gabor frames by exploiting knowledge about a related 
painless frame, to obtain "almost painless nonstationary Gabor frames". In Section [4] we 
provide examples based on the two introduced construction principles. 



2. Notation and Preliminaries 

Given a non-zero function g G -^ 2 (M), let gk,i(t) = M b iT ak g(t) := e 2mblt g(t — ak). M hl is a 
modulation operator, or frequency shift, and T ak is a time shift. 

The set Q(g, a, b) = {g^/ : k, I G Z} is called a Gabor system for any real, positive a, b. 
Q(g, a, b) is a Gabor frame for L 2 (M), if there exist frame bounds < A < B < oo such that 
for every / G L 2 (IR) we have 

A\\f\\l<J2\(f>^i)\ 2 <B\\f\\l (1) 

To every Gabor system, we associate the analysis operator C g given by (C g f)k,i = (f, gk,i), 
and the synthesis operator U 1 = C*, given by U y c = J2k iez c k,ilk,i for c G £ 2 . The operator 
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S 9il associated to Q(g, a, b) and Q(j, a', b 1 ), where S g>1 = U^Cg reads 

S 9,yf = ^2(f,9k,ihk,i- 

k,i<=z 

The inequality (pj) is equivalent to the invertibility and boundedness of the frame operator 
S g ,g of G(g,a,b). 

The analysis operator is the sampled short-time Fourier transform (STFT). For a fixed 
window g G L 2 (R), the STFT of f £ L 2 (R) is 

V g f(x,u>)= [ f(t)e- 2 ^ t gJt^)dt=(f,M UJ T x g). 
Jr 

Setting (x,co) = (ak,bl), leads to V g f(ak,bl) = (C g f)k,i- 

When working with irregular grids, we assume that the sampling points form a separated 
set: a set of sampling points {a^ : k G Z} is called (^-separated, if \a,k — a m \ > 5 for a*, a m , 
whenever k 7^ m. xi w iU denote the characteristic function of the interval /. 

A convenient class of window functions for time- frequency analysis on L 2 (IR) is the Wiener 
space. 

Definition 2.1. A function g G L°°(M) belongs to the Wiener space IV (L 00 ,!' 1 ) if 

\\g\\w{L^^) ■= 5^ ess su p*gqI^(* + k )\ < 00 ' Q= [M] ■ 

fcez 

For 5- G ^(L 00 , £ x ) and 5 > we have [15] 

esssup tgM ^|#(t-<5A;)| < (1 + (T 1 ) ||#|| • (2) 

kez 

In dealing with polynomially decaying windows, we will repeatedly use the following 
lemma. 

Lemma 2.2. For p > 1 the following estimates hold: 

(a) Let 5 > 0, then 

00 

5^(1 + sky p < (1 + 5)- p {5- 1 +p)(p- . 

k=l 

(b) Let {at : k E Z} <zR be a 5— separated set. Then 

esssu ta ^(1 + |t - a fe |)- p < 2 (l + (1 + ^(r 1 +p)(p - l)" 1 ) . 

fcez 

Proof. To show (a) we write 

00 00 00 „ 

5^(1 + = (1 + + ^(1 + = (1 + 5)- p + i 1 + Sk )' P dt ■ 

k=l k=2 k=2 J m+k 



3 



for t G [kj k + 1], we have St < 5(k + 1) which implies that 1 + 5(t — 1) < 1 + 5k. Therefore, 

oo « oo « /*oo 

V / (1 + 5A;)- p dt<J2 (1 + S(t- l))- p dt = (1 + 5(t - l))- p dt 

k=2 J[0,l]+k k=2 J[0,l]+k J2 

= (i + 5)-^ 1 r 1 ( P -i)- 1 , 

and the estimate follows. 

To prove (b), fix t G R. Since \a k — a\\ > 5 for k ^ I, each interval of length S contains at 
most one point t — a k , k G Z. Therefore we may write t — a k = Sn k + x k for unique n k G Z 
and x k G [0, 5), and by the choice of S, we have n k ^ rii for k ^ I. We assume, without loss 
of generality, that n k = for k = 0, and we find that 

]T(1 + |t - a k \)- p = ]T(i + \5n k + x k \)- p 
fcez fcez 

<1+ ^ + Sn k + x k )~ JJ + (^ + Sn k -x k )- p 

k£Z;n k >0 fceZ;n fe >0 

<1+ + J2 (l + tn k -5)- p 

k£Z;n k >0 fceZ;n fc >0 
oo oo 

< 1 + ^(1 + 6k)- p + ^(1 + S(k - l))- p 

k=l k=l 

= 2 f 1 + JO 1 + 5A; )~ P J < 2 (i + (i + 1 + P )(p - i)- 1 ) . 

The last expression is independent of t, and the claim follows. □ 
Remark 1. When the set A = {a k : k G Z} C E is relatively 5— separated, meaning 

rel(.A) := max #{A D ([0, 5] + t)} < oo , 



then the estimate (b) in Lemma 2.2 becomes 

esssu PteM ^(l + |t-a fe |)- p <2rel(^) (l + (1 + <5)- p (cT 1 + p)(p - l)" 1 ) . 

fcez 

Notice, that for a separated set, rel(^4) = 1. 

3. Nonstationary Gabor frames 

Nonstationary Gabor systems provide a generalization of the classical Gabor systems of 
time-frequency-shifted versions of a single window function. 

Definition 3.1. Let g = {g k G L 2 (IR) : k G Z} be a set of window functions and let b = 
{b k : k G Z} be a corresponding sequence of frequency-shift parameters. Set g*:,z = Mb k ig k . 
Then, the set 

Q(g,h) = {g k y. k,leZ} 
is called a nonstationary Gabor system. 
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Note that, conceptually, we assume that the windows gk are centered at points {a^ : 
k G Z}, in direct generalization of the regular case, where guit) = g(t — ak) for some time- 
shift parameter a. In this sense, we have a two-fold generalization: the sampling points 
can be irregular and the windows can change for every sampling point. We are interested 
in conditions under which a nonstationary Gabor system forms a frame. We first recall 
the case of nonstationary Gabor frames with compactly supported windows, see pj and 



http : //www . univie . ac . at/nonstatgab/ for further information. 



3.1. Compactly supported windows: the painless case 

Based on the support length of the windows g k , we can easily determine frequency- shifts 
parameters bk, for which we obtain a frame. The following result is the nonstationary version 
of the result given in [7]. 

Proposition 3.2 Q2]). Let g = {gk G L 2 (IR) : k G Z} be a collection of compactly supported 
functions with \supp g%] < 1/bk- Then Q(g,h) is a frame for L 2 (IR) if there exist constants 
A > and B < oo such that 

A<G (t) = J2 b k 1 \9k(t)\ 2 < B a.e. . 

k€Z 

The dual atoms are then 7fc,;(t) = M^Gq 1 (t)gi e (t). 

Remark 2. An analogous theorem holds for bandlimited functions gk- 

The above theorem follows from the fact that the frame operator associated to the col- 
lection of atoms described in the theorem can be written as 



S 9 J(t) = J2 b k 1 \9^)\ 2 f(t) a.e. 



fcez 



The diagonality of the frame operator in the painless case is derived from a generalized 
Walnut representation for the frame operator S* 9j9 of nonstationary Gabor frames. In the 
next section we will see that this representation immediately implies diagonality of S 9} g under 



the assumptions of Proposition 3.2 



3.2. A Walnut representation for nonstationary Gabor Frames 

Let us now consider nonstationary Gabor systems Q(g, b) and £(7, b), with all windows 
gk and 7^ in W^L 00 ,? 1 ). The operator associated to (?(g,b) and £7(7, b) reads 

W= Y^^ M ^k)M lbhlk . (3) 



Proposition 3.3. The operator S gn in Q admits a Walnut representat 



ion 



s 9 ,jf = E G l'J ■ > where g m W = b ^ - l K l )ik{t) , (4) 
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for f G L 2 (IR). Moreover, its operator norm can be bounded by 

( \ x l 2 ( \ 1//2 

\(S g ,yf,h)\ < (^sup(l + 6 A : 1 )||7 fe || H / (L ^ / i ) J (esssup tgR 2j# fc (t)|J 

k&L k&Z 

sup(l + ^ 1 )||^ fc || w(i o ^i ) ) ( ess sup igK /\lkif) I ) H/H2INI2 (5) 
key. / V _ — _ / 



for all f,h£ 1/ 

Proof. First assume that f,h e L 2 (IR) are compactly supported. Since (f,Mib k gk) 
(fgk)(lbk) we can write 5^ as 

S fll7 /(*) = E (7S)(i6fc)M % 7 fc (t) = ^m fc (t) 7fc (t) 



(6) 



where m fc (t) = J2i eZ (fgk)(lb k )e 2mlbkt , for every 6 Z. The functions are 6 fc 1 periodic 
and by the Poisson formula can be written as 



m k 



(t) = b k 1 J2(f9k)(t-lb k 1 ). 



(7) 



Therefore, substituting d7]) in (JHJ yields the Walnut representation. 

We next prove the boundedness ([5]). In the following chain of inequalities, we will use 
Cauchy-Schwartz inequality for sums and integrals and, since all summands have absolute 
value, Fubini's theorem to justify changing the order of summation and integral. We thus 
find 



\(S„„f,h)\ = |( E K 1 9 k (--lb k 1 h k (-)f(--lb k 1 ),h 

< E 6 ^ 1 / Mt-ib^umnt-iK^imidt 
E 6 * 1 IfMt-iK^uimfit-ib^dt] 1 \\\g k {t-ibt)\\ik{t)\\h{t)\ 2 dt 



k,iez 



< 



E 6 ^ / \9k(t)\hk(t + lb^)\\f(t)\ 2 dt 



imi^^Mmikit-ib-^idt 



kl&L 



1/2 



1/2 



E^/It^H^-^ 1 )!!^)! 2 ^ 
\Kt)\ 2 Y, h k l hk(t)\\gk(t-ib k 1 )\dt 



1/2 



1 1/2 



kl&L 



(8) 
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The first term in the last expression can be bounded as follows 



\f(t)\ 2 £ K'MVUit - ib?)\dt = f £|t*(* - ib k l )\\f{t)\ 2 \g k (t)\dt 

<J2 (^esssup^E^C*-^ 1 )!] / \f(t)\ 2 \g k (t)\dt 
<sup (ft^esssup^El^C*-^ 1 )! J £ / l/(*)| 2 |<7k(*)l<** 
< sup b^ 1 ess sup tm ^2\-f k (t - lb^)\ esssup teM £|# fe (t)| 



11/112- (9) 



Using relation (J2]) for ess sup igK Xwezl7fc(* — ^ k l )\ anc ^ substituting (|9| into ([8| yields ([5]). 
The estimate for the second term in Q is obtained analogously. By the density of compactly 
supported functions in L 2 (IR), the estimate holds for all of L 2 (IR). □ 



Remark 3. From (|8j) in the proof of Proposition 3.3, it follows that the operator S ga is also 
bounded by 

\{S g;J f,h}\ < ( esssup teIR £ K^kit - Ib^lhki^l 

k,lez 

esssup teR £ ^iTfc^-^fc 1 )!!^*)! 



2 ' 4 2 



Remark 4. In the case of a frame operator S^g, the Walnut representation Q becomes 



S 9 J(t) = £ 6 ^(* " ^QkWit - lb~ k l ) a.e. , (10) 



fc,zez 

and the above bounds reduce to 



K^W^I < (sup(l + 6 fc 1 )||pfc||iv(L«»^i)) fesssup teK £|£ fc (t)|)||.f|| 2 ||/i|| 2 
|<<W^)I< (esssup^E^I^^-^II^Wl) ll/IWMa 



Remark 5. Note that in the painless case of Theorem |3.2[ the frame operator reduces to the 
multiplication operator Sg )9 f = J2kez^ko ' f = ^o • /■ 

Remark 6. In the standard setting of Gabor frames, i.e. gk{t) = g(t — ak) for fixed a > 0, 
and b k = b for all fceZ, the above bound reduces to the well know bound 

\(S 9 , 9 f,h)\< (l + a-^Cl + r^ll^ll^^ll/llall/illa 
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3.3. Existence of nonstationary Gabor frames 

For windows that are neither compactly supported nor bandlimited, we are interested 
in the existence of frames of the form Q(g, b) and in the construction of the involved param- 
eters. The following theorem derives a sufficient condition for the existence of nonstationary 
Gabor frames and shows that this condition can be satisfied. 

In this and the subsequent sections, [bi,bu], \pi,Pu], [Cl,Cu] are compact intervals of 
positive real numbers. 

Theorem 3.4. Let g = {g k G W(L°°, I 1 ) : k G Z} be a set of windows such that 

i) for some positive constants A ,B 

0<A < J2\9k(t)\ 2 <B <oo a.e. ; (11) 

ii) for all k G Z, the windows decay polynomially around a 5 -separated set {a k : k G Z} 
of time-sampling points a k 

\9k(t)\ <C k (l + \t-a k \)- p \ (12) 
where p k G [p L ,Pu] CR, p L > 2 and C k G [C L , Cu}. 

Then there exists a sequence {b k } ke z, such that for b k < b k , k G Z, the nonstationary Gabor 
system G(g, b) forms a frame for L 2 (IR). 



Proof. Let / G L 2 (IR). Applying (10), we write 



w/>= / Y, b k 1 \9^ t )\ 2 \f(t)\ 2 dt+ f e J2 b k 1 9dt)g k (t-ib k i )f(t-ib k i )fxt)dt 

^ R fcez ^ R zez\{o} fcez 

Using similar arguments as in the derivation of (J8l), we obtain 

/ E Y, b ?9k(t)9k(t - ib k i )f(t - ib^wjdt < 

<lesssup teR E E 6 * 1 ^*^!!^ - ' 6 * 1 
\ zez\{o} fcez 

< ^{^'} E ( esssn Ptm^2\9k(t)\\g k (t-lb 



■fc 1 )! 11/111 



R 



Therefore, lower and upper frame bounds are obtained from 

(S g J,f) ||/|| 2 2 >min{6^} fessinf teR El^WI 2 " m ^ fcez{ g% ) (13) 
k& V k& mm k& {b k } ) 

(Sg, g f,f) 11/112 2 ^maxjfefc 1 } esssup^ElffcWI 2 + R 
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We need to construct a sequence of bk, fceZ, such that for all / G £ 2 (M), ( |13| ) is bounded 
away from zero. 

Let e < Cl and consider the sequence of frequency shifts b k = (ff-) 1 ■ Then min fcgZ {&k > 
(Clc- 1 ) 1 ^, max^^T 1 } < (C^e" 1 ) 1 ^ and 

maxfcgzj^ 1 } < c x/ PL( j-i/pv e ifp v -xfPL (x4) 
min fceZ {6^ 1 } - ^ L 



Since (1 + |a; + y|) p < (1 + |£|) p (l + \y\) p for x, y G R and p > 0, using (12), we have, for 
some /i with p^ — 2 > /i > 0: 

- < C 2 k (l + \t- o fc |)-»(l + |t - a, - /6-i|)-»+(i+") 

< C 2 (l + |t - a fc |) (1+ ^(l + |/|& fc T Pfc+{1+M) 

< c 2 (l + |t - a k \)-^ 6P*-d+") 

= cl +{1+ ^ /Pk (l + \t- a k \)-^ \l\-Pk+(i+^) e i-(i+M)M . 

Hence, 

fcez fcez 

< maxE k \1\-pl+(i+^ e i-(i+M)/Pt y^(! + u _ qJ)-^) 



< maxE fc e i-d+f)/Pt 2 (1 + (1 + + 1 + ^p" 1 ) , (15) 



where the last estimate f ollow s from Lemma 2.2 (b). Summing the expression (15) over 



I G Z\ {0} using Lemma 2.2 j , we see that R, as a function of e, behaves like e 1 TR-aOM^ 



i.e. i? ~ e 1 C^w/p^, and i? tends to for e — >• 0. Moreover 



maxfcgzj&fc 1 } 
minfcgzjfe^ 1 } 



i? ftj e 1 -( 2 +M)/PL + VpL 



can be made arbitrarily small by choosing e small since 1 — (2 + /i)/pL + 1 /W > 0. Therefore, 
if e is such that for b° k := (^) 1/pk , ^^^j R < A), then {M lbk g k } Kl& is a frame for all 

h<b° k . " & " □ 

For completeness, we state the equivalent result for analysis windows g k with polynomial 
decay in the frequency domain. 

Corollary 3.5. Let g = {g k G L 2 (IR) : ^ G W-^L 00 ,^ 1 ), k £ Z} be a set of windows such 
that 

i) for some positive constants A ,B Q 

0<A < ^|^(t)| 2 < B < oo a.e. ; (16) 
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ii) for all fcgZ, the windows decay polynomially around a 5 -separated set {b k '■ k G Z} 
of frequency-sampling points b k : 



\g k (t)\ < c k (i + \t - b k \y 



-pk 



(17) 



where p k and C k are chosen as in Theorem 3.4 



Then there exists a sequence {a k } k ^z, such that for a k < aP k , k G Z ; the nonstationary Gabor 
system {Ti ak g k : forms a frame for L 2 (IR). 



3.3.1. Nonstationary Gabor frames on modulation spaces 

Modulation spaces, cf. [T3], US], are considered as the appropriate function spaces for 
time-frequency analysis and in particular, for the study of Gabor frames. By their definition, 
modulation spaces require decay in both time and frequency. Under additional assumptions 
on the windows g k , the collection Q(g, b) is a frame for all modulation spaces M p , 1 < p < oo. 

Proposition 3.6. Let Q(g,b) be a frame for L 2 (K) satisfying the uniform estimate 

\V*g k (x, u)\ < C(l + \x- a fc |)- r " 2 (l + \u\)- r ~ 2 , r > 2 (18) 

where cf) is a Gaussian window. Then the frame operator S is invertible simultaneously on 
all modulation spaces M p for 1 < p < oo. 

Proof. Notice, that 

\V4,g h j,{x,u)\ <C(l + \(x- a k ,uj- lb k )\y r - 2 , 

since (1+ \x\ + |w|) _r > (1 + |x|)~ r (l+ |u;|)~ r and the weights (1 + \(x, uj)\) r and (1 + \x\ + |o;|) r 
are equivalent. 



A result on Gabor molecules [16] states that, if an L 2 -frame {g z : z — {z\, z%) G Z C ]R 2 }, 
where Z is separable, satisfies the uniform estimate \y^g z {x,uS)\ < C(l+\(x— zi,u>— z 2 )\)~ r ~ 2 , 
then the frame operator Sf = Ylzezifi 9z)9z is invertible simultaneously on all M p for each 



1 < p < oo. The result hence follows from condition (18). □ 
3.4- Constructing nonstationary Gabor frames 



Theorem 3^ shows that for windows with sufficient uniform decay, nonstationary Gabor 
frames can always be constructed by choosing sufficient density in the frequency samples. 
In the present section we assume the existence of a certain nonstationary Gabor frame and 
explicitly construct a new frame by exploiting the prior information about the original one. 
This is a situation of practical relevance, since we may often be interested in using windows 
that decay fast and are negligible outside a support of interest. In particular, we will use 
the fact that painless nonstationary Gabor frames are easily constructed and deduce new 
frames from painless frames. The new frames thus obtained will be called almost painless 
nonstationary Gabor frames. 

We will subsequently assume b k e [£>£,, bu] for frequency-shift parameters and we let C k E 
[Cl,Cu] and p k G [pl,Pu] with pi > 1 for the constants involved in the decay assumptions 
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for the analysis windows. We then work with the following constants that depend on the 
separation of the sampling points, the decay of the windows and the frequency-sampling 
parameters: 

7? 1 , ^+Pl , „ -. bu+pu 

El = 1 + — rr -, — and E 2 = 1 + — rr ; r (19) 

(i + «j)^(p L -i) (l + r^^-i) 1 ; 

We first consider nonstationary Gabor frames obtained by perturbation of a known frame. 
This result is in the spirit of similar results for regular Gabor frames [3j E]. 

Proposition 3.7. Let {a k : k G Z} be a S-separated set and Q(h, b) a nonstationary Gabor 
frame with frame bounds Ah, Bh and frame operator Sh,h- Let G L 2 (IR) be a set of windows 
such that for all k G Z and /or almost all t 6l 

|</*(t) - &*(t)| < C fe (l + |t - a fc |)~ Pfe . (20) 

IfCu<y/A^for 

A = 46,; 1 • £ a • £ 2 , (21) 



£/ien £(g, b) a frame for L 2 (IR) with frame bounds A = Ah(l — \l C^XA h 1 ) 2 and B 



B h (l + y/C^XB^f. 

Proof. By applying Cauchy-Schwartz inequality, it is easy to see that, for Q(g,b) to be a 
frame for L 2 (IR), it suffices to show that ^2 k i GZ \{f, gk,i — h k ,i)\ 2 < f° r some R < A h , 

also cf. [6], Proposition 4.1.]. Then, frame bounds of Q(g, b) can be taken as Ah(l — a/ R/Ah) 2 
and B h (l + ^JR/B h ) 2 . 

To obtain the required error bound, we let ipk(t) '■— g k — h k and use the estimate given 
in (|9l): 



£l</> VOI 2 < sup (V esssu ta ^ fc (t - lb~ k l )\) (ess sup t6R |) ||/||£ . (22) 



The first term of the above in the above inequality is bounded by assumption (20) and 



Lemma |272| (b) : 

< Q/J> + \t- a k\r pk < Cc/^(1 + \t-a k \)-n<2CuE 



l • 



To bound the second term, note that ^ igZ (l + \t — a k — Ib^l) Pk is b k x — periodic, therefore 
we can simplify 



esssup^^l^^-/^ 1 )! < C fc ess sup t6[0)6 -i] J^(l + \t - Ib^l)' 

lei lei 
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-Pk 



Hence, by Lemma 2.2 we obtain for t G [0, b k 1 ]: 

oo oo 

E( x + i* - ^'D^ < i + Ea + 1* - lb k 1 \y pk + E^ 1 + * + "t 1 ) 

zez i=i z=i 

(1 + (i - 1)6^)-** 

oo oo 

= 2 £(l + i6-i)-» = 2(l + £(l + i6^)-») < 2£ 2 . 
Gathering all the estimates, we obtain 

EK/^^l^^r^i-^n/ii^o 



fc.i 



By assumption CV}A < A h , and the proof is complete. □ 



Using Proposition 3/T we next construct a special class of nonstationary Gabor frames by 
relying on knowledge of a painless nonstationary Gabor frame. We construct new windows 
which are no more compactly supported, but coincide with the known, compact windows 
on their support. We call the resulting new systems almost painless nonstationary Gabor 
frames. 

Corollary 3.8. Let g = {g k G W^L 00 ,^ 1 ) : k G Z} be a set of windows, and let I k be 
the intervals I k = [a* — (2&fc) -1 , a* + (26fc) -1 ] where {a k : k G Z} forms a 6-separated set. 
Assume that Q(h, b), where h k = g k Xi k , is a Gabor frame with lower frame bound Ah, and 
that for ip k = g k — h k , all k G Z and almost all tiriM. 

lN |-Pfc + n , -L (Oh, \- 1 ■ 



C k (l+t-a k - (2b k )- l )- p \ t>a k + (2b 



h. 



\4> k (t)\ < { 0> tel k ; (23) 

C fc (l - t + a k - (2b k )- 1 )-»> , t < a k - (2b k y l . 

If Cu < V A h X~ l for A = Ab] 2 ■ d^ 1 ■ E\ ■ E 2 , then Q (g, b) forms a nonstationary Gabor frame 
forL 2 (R). 



Proof. The proof follows the steps of the proof of Proposition |3.7| with small changes on how 
to approximate the terms in (22). First observe, that for any t G R and all k, we have 

\Mt)\ < C k [(1 + \t-a k - (26 fc )- 1 |)- pfc + (1 + \t - a k + (2b k )- 1 \)^] 

Since the frequency shifts b k are taken from a finite interval and the set {a k : k G Z} is 
5— separated, the sets T + = {a k + {2b k )~ 1 : k G Z} and T~ = {a k — (2b k )~ l : k G Z} 
are relatively 5— separated with rel(r) = rel(r + ) = rel(r~) = L(2&l5) _1 J. Therefore, by 
Remark [TJ it follows that 



Ei^wi ^ °u E^ 1 + i* - ak - ^y l \y pL + °u E^ 1 + 1* - % + ^r'lr 

fcez fcez 
< 4C c/ rel(r) (! + (! + 5)- pL (5~ l +p L )(p z - l)" 1 ) . 



Pi 
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Now, the expression XwezlV'fc^ ~~ ^fe 1 )! * s 1 — periodic. Let t G I k , then, by (23) 



J2(l + a k -(2b k )- 1 -t + lb k 1 ] 

_l>0 

+ J2^-ak-(2b k y 1 + t-lb k 1 )-^ 

oo 

< C k 2 ^(1 + lb k Y Pk < 2C fc (l + (1 + 6* 1 )- p *(6* + P*)(p* - i)" 1 ) , 
z=o 

where the last estimate follows from Lemma 2.2[ a). □ 



4. Examples 

We illustrate our theory with two examples. In both examples, we consider a basic 
window and dilations by 2 and |, respectively. Since the dilation parameters take only three 
different values, there are three kinds of windows, with support size 1/2, 1 and 2, respectively. 
Note that, while theoretically possible, sudden changes in the shape and width of adjacent 
windows turn out to be undesirable for applications, hence we only allow for stepwise change 
in dilation parameters. 

In the first example we consider a frame that arises as perturbation of a painless non- 
stationary Gabor frame. The perturbation consists in the application of a bandpass filter in 
order to obtain windows with compact support in the frequency domain. 

Example 4.1. Let h be a Hann or raised cosine window, i.e. h(t) = 0.5 + 0.5cos(27rt) for 
t G [—1/2,1/2], and zero otherwise. We construct a painless nonstationary Gabor frame 
by dilating h by 2 or |, respectively: Let {s k } k& i be a sequence with values from the set 
{—1,0,1} with the restriction that \s k — Sfc_i| G {0,1} to avoid sudden changes between 
adjacent windows. We then define corresponding shift-parameters by setting ao = and 



fl/t+i — «fc + 2 

a k+l = a k + 2 



-Sfc + 1 



3 
5 



if s k > s k+ i 
if s k = s k+ i . 



a k+1 = a k + 2 Sk+1 ■ - if s k < s k+1 . 
o 

The points a k , k G Z, form a separated set with minimum separation 5 = 1/3. Setting 
b k = 2 Sk and h k (t) = T ak \/2^h(2 Sk t) ) the system {M ibk h k : k,l G Z} forms a painless 



nonstationary Gabor frame with lower frame bound 
Let Q = 0.02 and be a bandlimited filter given by 



0.5. 



0.5 + 0.5 cos(27rrr ^ on its support [-Q/2, fl/2) 



We build new windows g k by convolving with h k . The windows g k '■= <p * h k are no more 
compactly supported. Since \4>(t)\ < Q(l + |t|)~ 3 , we rely on [201 Theorem 9.9] to deduce 
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the following bound, with C = \\hk\\oo%, h = [a* — 2 s ' k ,a k + 2 s *] and s' fc = + 1: 



(1 + (t - a fc ) - 2-<)- 2 - (1 + (t - a fc ) + 2- s ' fc )- 2 t > a k + 2-4 

2 _ (1 + (f _ a k ) + 2"4)-2 _ (i _ (f _ afe ) + 2-t)- 2 t g J fc 

(1 - (t - a k ) - 2-4)- 2 - (1 - (t - a fc ) + 2-4)-2 t < a fc - 2~< 



We obtain the joint bound, — h k {t)\ < Cu(l + \t — a k \) 2 by setting C k = C ■ (1 + 2 V 
for all k G Z and = max fceZ C fc = 0.0282 < y/AhX- 1 = 0.0768, with A as defined in pT 



Thus, by Proposition 3.7 {Mib k g k : k, I G Z} is a Gabor frame with a lower frame bound 
A = 0.2. 



Remark 7. Note that the construction presented in the Example |4.1| is of particular interest 
for constructing frequency- adaptive frames with windows that are compactly supported in 
time. This is a situation of considerable interest in applications, since it allows for real-time 
implementation with finite impulse response filters, cp. |llj . 

In the second example we construct a nonstationary Gabor frame by applying Corol- 



lary |3.8[ The windows of the new frame coincide with the windows of a painless frame on 
their support. The windows in this example are constructed in analogy to the windows used 
in scale frames, introduced in [2] to automatically improve the resolution of transients in 
audio signals. 

Example 4.2. As in the previous example, let s k G {—1, 0, 1} with \s k — s^-il G {0, 1} for all 
fceZ. We consider a sequence of windows g k that are translated and dilated versions of the 
Gaussian window g(t) = e ~ n(2M)2 \ g k (t) = T ak V2^g(2 s H) with a = and 





= a k + 2 


Sfe+i-l 


if 


Sk — S k+ i 








if 




a k+l 


= a k + l 


2-sk+i 


Sk > S k+ i 








if 






= a k + l 


2 Sk 


Sk < S k+ i 



Here, the {a k : k G Z} are separated with minimum distance 5 = 1/4. We arrange the 
windows as follows: after each change of window size, no change is allowed in the next step; 
in other words, each window has at least one neighbor of the same size. An example of the 
arrangement is shown in Figure [T} 

Let I k = [a k — 2~ Sk ~ 1 ,a k + 2~ Sk ~ 1 } and define a new set of windows by h k (t) = g k (t)xi k - 
Then {Mit, k h k : k,l G Z} with b k = 2 Sk is a painless nonstationary Gabor frame. By 
numerical calculations, its lower frame bound is = 0.1609 and ip k (t) = g k (t) — h k (t) can 
be bounded by 

v / 2^^(l/2)(l +t-a k - 2- Sk - 1 y 19 t> a k + 2~ Sfc " 1 

\Mt)\ < { o tei k 

s/2^g(l/2)(l + a k - 2~ Sk - 1 - ty 19 t < a k - 2~ Sfc ~ 1 . 



From Proposition 3.8 it follows that 4(56|)- 1 • ■ E 1 ■ E 2 = 0.0071 < A h , and {M lbk T ak g k : 



k,l G Z} is a nonstationary Gabor frame with a lower frame bound A = 0.1538. 
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Example of window arrangement 
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Schematic depiction of resulting time-frequency tiling 
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Figure 1: An example for the arrangement of dilated windows in Example 2 
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